Persistent spin currents in mesoscopic Haldane-gap spin rings 
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Using a modified spin-wave approach, we show that in the presence of an inhomogeneous magnetic 
field or an in-plane inhomogeneous electric field a mesoscopic antiferromagnetic Heisenberg ring 
with integer spin (i.e., a Haldane-gap system) exhibits a persistent circulating spin current. Due to 
quantum fluctuations the current has a finite limit on the order of {—g^B)c/ L at zero temperature, 
provided the staggered correlation length ^ exceeds the circumference L of the ring, in close analogy 
to ballistic charge currents in mesoscopic normal-metal rings. Here c is the spin-wave velocity, g 
is the gyromagnetic ratio, and hb is the Bohr magneton. For ^ <^ L the current is exponentially 
suppressed. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.30.Ds, 73.23.Ra 



I. INTRODUCTION 

In a recent paper, ^ henceforth referred to as I, we 
have shown that the mesoscopic persistent current in a 
normal-metal ring threaded by an Aharonov-Bohm flux 
has an analog in spin rings: specifically, we have calcu- 
lated the persistent spin current Is in thermal equilibrium 
that circulates a ferromagnetic Heisenberg ring subject 
to a crown-shaped magnetic field with magnitude 
Within linear spin-wave theory we have shown that the 
associated magnetization current /m = {g^BlK)Is (where 
g is the gyromagnetic ratio and /xb is the Bohr magneton) 
is of the form.^ 
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FIG. 1: Classical ground state of an antiferromagnetic Heisen- 
berg ring in a crown-shaped magnetic field; mf and mf are 
the directions of the spins on sublattices A and B. 



L has a sawtooth shape as a function of the geometric 
flux 



where L is the circumference of the ring, Ck is the energy 
dispersion, and Vk = h^^dek/dk is the velocity of the 
magnons, \h\ = giJ,B\B\, and T is the temperature. In 
the presence of a crown-shaped inhomogeneous magnetic 
fleld the direction of the magnetization in the classical 
ground state covers a flnite solid angle Qm on the unit 
sphere as we move once around the ring. Then the Bloch 
wave vectors are quantized as k„ = ^(n + ^^), where 
in the continuum approximation the allowed values of n 
are n = 0, ±1, ±2, . . .. 

The magnetization current vanishes for T — > 0. 
Physically, this is due to the fact that no quantum fluc- 
tuations are present in a ferromagnet, so that at T = 
there are no magnons to carry the spin current. It is thus 
tempting to speculate that the analogous current for an 
antiferromagnetic Heisenberg ring will be flnite even at 
T = due to quantum fluctuations. In this work we shall 
show that this is indeed the case and present a quantita- 
tive calculation of the current in Haldane-gap antiferro- 
magnets (i.e., antiferromagnetic Heisenberg rings with in- 
teger spin S) using a modifled spin-wave theoryi^*^ Our 
main result is that the ground state of a Haldane-gap spin 
ring subject to an inhomogeneous magnetic fleld supports 
a flnite magnetization current, which in the limit where 
the staggered correlation length ^ is large compared with 
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where c is the spin-wave velocity. Here f2 is the solid angle 
on the unit sphere traced out be the local Neel vector as 
one moves once around the ring. Such a state can be pro- 
duced by an inhomogeneous magnetic fleld, as depicted in 
Fig. n Similar to the case of a ferromagnet discussed in 
I, the magnetization current is carried by magnons which 
are subject to mesoscopic interference due to the geomet- 
ric phase associated with the inhomogeneous nature of 
the classical ground state. Due to quantum fluctuations, 
the ring is endowed with an electric dipole moment even 
in the ground state. The absence of true long-range order 
in one-dimensional Heisenberg antiferromagnets leads for 
integer spin 5 to a finite spin-correlation length ^ and a 
Haldane gap of the order of fic/£^ between the ground- 
state energy and the lowest triplet excitation. These 
features are correctly captured by a modified spin-wave 
theory'^'^ which we use here to calculate the spin cur- 
rent in Haldane-gap systems. Note that this approach 
is only appropriate for integer S, where the low-energy 
excitations can be viewed as renormalized spin waves. 
In contrast, for half integer 5* the spectrum is gapless 
and spin correlations decay algebraicallys^ The elemen- 
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tary excitations are then spinons, so that the modified 
spin-wave theory does not correctly reproduce the low- 
energy physics. In this case the effective low-energy the- 
ory is a Tomonaga-Luttinger model. Recently Meier and 
Loss'' used such a model to discuss the spin current in 
5=1/2 antifcrromagnetic spin chains for a two-terminal 
geometry. 

While Eq. ^ is the bosonic analog of the persistent 
charge currentSiS in a normal one-dimensional metal ring 
threaded by an Aharonov-Bohm flux 0, Eq. J^J is for- 
mally identical with the persistent charge current in a 
ballistic metal ring at zero temperature. Recall that for 
spinless fermions at constant chemical potential /x the 
charge current in the ballistic regime can be written as^ 



generalize the spin- wave approach, we then apply a mod- 
ified spin-wave theory, where the absence of long-range 
order in a one-dimensional antiferromagnet is taken into 
account. In Sec. IIIII an antiferromagnet in an inhomo- 
geneous electric field is considered. In Appendix A, the 
classical ground state is obtained explicitly for the sim- 
ple geometry of a crown-shaped field, and Appendix B 
contains some mathematical details of the calculation of 
the magnetization current. 

II. SPIN- WAVE THEORY 
A. Spectrum 



k 

Here is the energy dispersion and Vk is the correspond- 
ing velocity of an electron in the state with wave number 
k. Due to the Aharonov-Bohm flux the wave vectors are 
quantized as fc„ — "^{n — (p/cj)^), where 0o = hc/e is the 
flux quantum. For free fermions at zero temperature Eq. 
P|) reduces to a sawtooth function. For an even number 
of electrons at T = one obtains 

O-^)' ^<i<- 

where vp is the Fermi velocity. With the replace- 
ments e 5/iB, Vp — > c, and 0/0o r2/27r the 
zero-temperature charge current (0J is formally identi- 
cal with the zero-temperature magnetization current (0) 
in a Haldane-gap spin ring. Finite temperature, disorder, 
and phase-breaking scattering all have a similar effect on 
the persistent charge current in one dimension. For a 
weak perturbation, they smooth the discontinuity around 
(j) = and with increasing strength higher harmonics are 
exponentially suppressed such that a sinusoidal shape is 
approached. In the limit of a very strong perturbation 
the current is exponentially suppressed with the relevant 
length or energy scale, i.e., the current becomes propor- 
tional to exp(— T/T*), exp(— L/L^), or exp(— L/L^) un- 
der the influence of a nonzero temperature, strong dis- 
order or strong inelastic scattering, respectively. Here 
T* is the temperature scale associated with the discrete 
level spacing, is the localization length, and is the 
phase-coherence length. Below we show that in the case 
of spin currents in Haldane gap systems the correlation 
length ^ plays a similar role: for ^ <C i the magneti- 
zation current is exponentially suppressed and becomes 
sinusoidal. 

The rest of the paper is organized as follows. In SecHll 
we apply linear spin-wave theory to an antiferromagnet 
subject to an inhomogeneous magnetic field. After cal- 
culating the spectrum of magnon excitations, we obtain 
the spin current as a gauge-invariant derivative of the 
flux-dependent part of the free energy and discuss the re- 
sult for T = in naive spin-wave theory. To justify and 



We start with a nearest-neighbor antiferromagnetic 
(J > 0) Heisenberg ring in an inhomogeneous magnetic 
field 

N 

where Si are localized spin operators at sites h of a one- 
dimensional lattice with lattice spacing a = L/N, and 
Sf = S{S+1). Periodic boundary conditions Sj\r+i = Si 
are used, hi — gfi^Bi, and is assumed to be even. For 
a spin- wave expansion the classical ground state has to be 
determined by minimizing the energy in Eq. (jSJ with Si 
replaced by Srhi, where rhi is a unit vector. To consider 
fluctuations, the spin operators are then represented by 
bosonic creation and annihilation operators in the stan- 
dard way, 

S^-m, = S- bib, , S, ■ e+ = V2Sb,[l + 0{S'')] , (6) 

where = + ie^, and and ef are two unit vectors 
in the plane perpendicular to m,;, such that {e}, ef, m^} 
is a right-handed local basis in spin space. There is a 
local gauge freedom in the choice of the transverse basis 
vectors e} and ef which can be arbitrarily rotated around 
rrii. With the notation of I the spin- wave Hamiltonian 
then reads 

i 

- (1 - m, • m,+i)e*("— + h.c] 

- 2mi ■ mi+i{blbi + bl_^_^bi+i) + hi ■ mi6j6,;| , (7) 

where iOi^j is the angle of rotation around rrii that takes 
rrii X rrij into e}. The transformation to another local 
right-handed triad, associated with site i and the bond 
{ij) according to {ej, ef = rrii x 'frli]^ then reads 

e±=e±-.-.e,f . (8) 

Equation {Tj) is manifestly invariant under the U(l) gauge 
transformation LOi-^j — > uJi-tj -l- a^, bi e^°'^bi. We will 
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now assume a sufficiently large ring so that the direction 
of the magnetic field varies only slightly on the scale of 
the lattice spacing a. The classical ground state then 
locally resembles a Neel state and the local Neel vector 
rii = (— l)'+^mi varies smoothly as a function of posi- 
tion on the lattice and is oriented almost orthogonal to 
the local direction of the magnetic field. We thus have 
rrii ■ rrii+i = — 1 + 0{1/N) and the terms involving the 
combination in Eq. |(7)) can be neglected to leading 

order in In this approximation a local spin devia- 

tion follows the direction of the classical ground state as 
it moves around the ring. It picks up a geometrical phase 
leading to interference of its wave function in close anal- 
ogy to Aharonov-Bohm interference in charge transport. 
The effect of the inhomogeneous field can be incorporated 
via the gauge transformation 



(9) 



The new boundary conditions then become in the bosonic 
language 



bi+N = e'^^^bi , ft = ON+i 



(10) 



where the upper/lower sign is valid for sublattice A/B 
(odd/even i). The resulting quadratic bosonic Hamil- 
tonian is standard for an antiferromagnetic ring with 
nearest-neighbor interactions 

Hafm = ^ [ - JS{bA+i + blbl^,) + {2JS + hs)b\b.] , 

i 

(11) 

except for the twisted boundary condition H1U|) . An ad- 
ditional staggered field hs in the direction of the classical 
ground-state vectors has been introduced as a technical 
tool for the discussion in Sec. Ill CI Hafm is diagonalized 
as usual by first performing Fourier transformations with 
different signs on the two sublattices: 



a-k 



i£A ieB 
where the allowed wave vectors are given by 

kn^j^[n+-), n^O,...,--l. (13) 

The diagonal form of Hafm is then achieved by the Bo- 
goliubov transformation 



_ cosh^fe sinhfffe W "fc , . .^ 

b[ ) ^[ sinh^fc cosh(?fe j [pl I ' ^'^> 



with 



tanh(20fc) = ^ hs = hs/2JS . (15) 

1 + hs 



The diagonal Hamiltonian contains constant terms due 
to quantum fluctuations, 

Hafm = ^ e^ialak + filh + 1) - NJS{\ + hs) , (16) 



where the quasiparticle energies are given by 



2JS'A/A2-hsin^(A:a), ^^ = hs{hs + 2), (17) 



and the free energy is obtained from 



F^^ifl) = 2T^ln [2sinh 



2T 



NJS{l + hs). (18) 



Thus we have shown that to leading order in 1/A^ ther- 
modynamic quantities depend on the inhomogeneity of 
the field only via the single phase ^l. Geometrically fl is 
the anholonomy associated with the parallel transport of 
a vector orthogonal to the local Neel vector jii in anal- 
ogy to the ferromagnetic case. To see this more clearly, 
we consider for each bond (ij) two additional sets of lo- 
cal right-handed triads containing the Neel vector rii in- 
stead of rrii. These triads are given by {e^ = e^,e| = 
(— l)*+^e|, n^} and {e^,e| = m x nj,ni}, and are re- 
lated by a rotation around rii. For the associated spher- 
ical vectors this reads 



(19) 



where the rotation angles are given by 

=iTT + {-ly+^uj^^j for j = i±l. (20) 
We can now express ft as 



N 



i+l — l^i+l- 



mod 27r . 



(21) 



which is of the form obtained in I for the ferromagnet. fl 
is thus the anholonomy of a vector orthogonal to the lo- 
cal Neel vector that is transported around the ring by 
discrete rotations around rii x Tii+i- Alternatively, a 
continuous parallel transport can be used around a path 
of geodesies connecting the unit vectors rii on the unit 
sphere. O is therefore equal to the solid angle subtended 
by this closed path of geodesies. 



B. Magnetization current 

In I we have shown that the U(l) gauge symmetry as- 
sociated with the choice of the local transverse basis is 
connected with a conserved current that was identified 
with the longitudinal component of the spin current. For 
an antiferromagnet the longitudinal spin current is con- 
veniently defined in the direction of the local Neel vector 
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and can be written as a gauge-invariant derivative of the 
free energy 



= -(-1) 



i+1 



>i+i 



OF 

sw 



(22) 

where li^i+i = JSi x Si+i is the spin current from site 
i to i + Similar to the ferromagnetic case, the presence 
of a longitudinal spin current can be understood semi- 
classically as follows. The local spin deviation from the 
classical ground state is essentially perpendicular to 
and varies slightly from i to i + I. The spin deviations 
on neighboring sites are therefore in a plane that does 
not contain rii, so that their cross product appearing in 
7i^i+i has a nonvanishing component Is in the direction 
of rii. This spin current Is corresponds to a current of 
magnetic dipoles that are locally oriented in the direction 
of the Neel vector rii , which varies smoothly as we move 
along the ring. The spin current generates an electric 
dipole field which has the same form as discussed in I. 
For the magnetization current we obtain 



L 



E 



Cfe Pfe 



(23) 



where Ck = h~^dek/dk is the velocity of a magnon with 
wave vector k and nj, — l/[exp{ek/T) — 1] is the Bose oc- 
cupation factor. The extra factor of 2 is due to the two 
degenerate spin-wave modes. Equation H23|l is the an- 
tiferromagnetic spin analog of Eq. |j3Jl. It clearly shows 
that the magnetization current is carried by antiferro- 
magnetic magnons, which at this level of approximation 
are the only quasiparticles available for transport. The 
current has a finite limit, even for vanishing Bose occupa- 
tion factors, due to quantum fluctuations. From Eq. I|23|) 
the current is clearly seen to be a periodic function of fl, 
so that the finite momentum sum can be further evalu- 
ated using the Poisson summation formula. — Some de- 
tails of the calculation are given in Appendix For a 
vanishing staggered field in the zero-temperature limit we 
obtain the simple result announced in the Introduction, 



n 

1 — 

TT 



n 

< — < 1 

271 



(24) 



Here l';!^ — —gfiBc/L is the magnetization current car- 
ried by a single magnon with the spin-wave velocity 
c = Cj,_,o+ at the center of the Brillouin zone. The saw- 
tooth shape (see the solid line in Fig. O of the current 
in Eq. I|24(l is reminiscent of Eq. for charge trans- 
port. Indeed, for T = Eq. H23|) is formally equivalent 
to Eq. Q for charge transport when the Fermi edge is 
replaced by the lower edge of the magnon band. 



C. Modified spin- wave theory 

The usual spin- wave theory employed so far is inconsis- 
tent when zero modes appear. Although the spin current 
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FIG. 2: Magnetization current in a ring with 100 spins for 
different values of the energy-gap parameter A. The plots are 
produced by numerically evaluating Eq. I|23|l . For A = 10~^ 
the curve is indistinguishable from the approximate expres- 
sion in Eq. 1281 , and Eq. 12911 provides a good approximation 
for A = 2 X 10"^ 



remains finite, the sublattice magnetization diverges in 
the limit ^ 0. This failure is related to the absence of 
long-range order in one-dimensional systems. It can be 
resolved by a modified spin-wave theory which was first 
used by Takahashi for a one-dimensional ferromagnet^ 
and then extended to various spin systems without long- 
range order, including antiferromagnetsi^ The constraint 
that is introduced in these theories was recently shown 
to follow naturally from a calculation at constant order 
parameter^ In the present context, we introduce the ad- 
ditional constraint 



(25) 



which suppresses Neel order on average. This constraint 
is enforced via the staggered field hs in Eq. Ullfl which 
acts as a Lagrange multiplier. The expectation value in 
Eq. H25|) can be evaluated from dF^^jdhs, yielding the 
self-consistency condition 



TV ^ dhs 

k 



S- 



(26) 



Although the self-consistently determined hg is itself a 
periodic function of the geometric flux the leading 
order for large is a constant and can be determined 
by replacing the sum in Eq. H26|) by an integral. For 
T = the solution of Eq. (|26|l yields the Haldane gap 
2JSA, which is inversely proportional to the staggered 
correlation length ^ ^ 



A 



V2^ 



(27) 



The functional form of the magnetization current shows a 
crossover between the two qualitatively different regimes 
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^ ^ L and ^ ^ L (see Fig. In the former case 
A ^ 2tt/N and at most one wave vector can be in the 
region — A < A: < A where the dispersion relation de- 
viates strongly from the dispersion in the limit A = 0. 
When the contribution from this single wave vector is 
taken into account separately, we obtain 



7°" 



sm{2n/N) 



-T:<n<TT. (28) 



In the case A = this reduces to Eq. (|24ll . provided 
^ 1. Thus, the effect of a finite A is to remove the 
discontinuity at = 0, 27r. On the other hand, in the 
limit ^ <C i many k values are affected by the energy gap. 
An analytic result for the spin current can nevertheless 
be derived as described in Appendix IbI We obtain the 
scaling form 



L 



I'm V7^V^/2e 



1/2 



cxp 



L 

Tie 



sin(r2) , (29) 



implying that the sinusoidal magnetization current is ex- 
ponentially suppressed in the bulk limit, L ^ ^. 



III. ELECTRIC FIELD 

Moving magnetic dipoles represent an electric dipole 
moment^^ and are therefore affected by electric fields. 
Due to this relativistic effect, which is essentially equiva- 
lent to spin-orbit coupling, the magnetic moments ac- 
quire an Aharonov-Cashcr phasejAi For localized spin 
systems described by a Heisenberg Hamiltonian, the elec- 
tric field can be taken into account phenomenologically 
by a substitution in the interaction term. 



Si • Sj 



(30) 



as long as the electric field varies only weakly on the scale 
of the lattice spacing. Here, 
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^/ dlxE{l), 



(31) 



and e^^ denotes the SO (3) rotation matrix acting on a 
vector m according tc^'' 

e^^m = e{e-m) + {9xm)sme~9x{9xm)coae , (32) 

with = 96. For ferromagnetic coupling, inhomoge- 
neous electric fields can lead to persistent magnetization 
currentgi^ and a spin analog of the Hall effect was also 
shown to exist in electric fields i 

We now consider the antiferromagnctic ring in an elec- 
tric field in the x-y plane, e.g. produced by a charged 
line in the z direction (see Fig. 13). The rotation vectors 
— ^4,1+1 are then all parallel to the z axis and 
the Hamiltonian for vanishing magnetic field becomes 



(33) 




FIG. 3: Antiferromagnetic Heisenberg ring in the electric field 
produced by a line charge. 



The classical ground state is easily shown to be a dou- 
bly degenerate Neel state with rrii = ±e^. The spin- wave 
expansion is thus straightforward. If a gauge transforma- 
tion is used to eliminate the phase factors, we again ob- 
tain the standard bosonic Hamiltonian Hafm of Eq. Hll|) 
with the boundary condition H10|l . where is replaced by 
the total Aharonov-Casher phase 



dl ■ \e'- X E(V)\ . 



(34) 



The spin current then only has a z component which 
can be written as a gauge-invariant derivative of the free 
energy. 



' on 



AC 



(35) 



This leads again to Eq. H23|) with fl replaced by Hac and 
all the results derived in the previous sections are also 
applicable in this context. 

It is also interesting to note that the situation of a 
radial electric field with = ^ and an additional 

homogeneous magnetic field tilted with respect to the z 
axis can be formally mapped onto a crown-shaped mag- 
netic field alone via the transformation 



e « 



(36) 



It would therefore be interesting to further investigate 
the combined effect of arbitrary inhomogeneous magnetic 
and electric fields on the produced spin currents to find 
situations that could be realized more easily in the labo- 
ratory for a possible experimental detection of the effect. 



IV. SUMMARY AND OUTLOOK 

In the last two decades, a lot of theoretical work fo- 
cused on persistent equilibrium currents in mesoscopic 
normal-metal ringsii2i24i& In the 1990s, the experimen- 
tal difficulties for the detection of the currents were 
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overcome and an oscillating magnetization as a func- 
tion of the magnetic flux was observed under vari- 
ous conditionsii£'iS*i2iSfliSiiSS In the ballistic regime, the 
main features of the experiments^' can be understood 
within a simple model of free fermions,'^" but in the 
diffusive regime a generally accepted explanation of the 
experimentsii^ is still lacking. 

In this work, we have shown that persistent magnetiza- 
tion currents are present in antiferromagnetic Heisenberg 
rings in inhomogeneous magnetic fields as well as in a ra- 
dial electric field. Quantum fluctuations lead to ground- 
state currents, but fluctuations in low dimensions also 
produce exponential damping when the circumference L 
of the ring becomes larger than the staggered correlation 
length. We have obtained explicit expressions for the 
current at T = in the two limits L ^ ^ and L 3> ^ 
within a modified spin-wave approach valid for integer 
spins, i.e., for Haldane-gap systems. The determination 
of the current for half-integer spin rings remains an in- 
teresting open problem. Also, we have only considered 
clean systems in this work, i.e., we have focused on the 
ballistic regime. Since for persistent charge current dis- 
order is known to be important, it would also be very 
interesting to consider persistent magnetization currents 
in the diffusive regime of disordered magnets. 

In the past few years a new field of research ( "spintron- 
ics") has emerged, where the spin degree of freedom is 
used as a medium to transport information^ For techni- 
cal applications (such as quantum computation) it is im- 
portant to control the quantum coherence in mesoscopic 
spin systems. The persistent spin current discussed in 
this work can be viewed as a quantitative measure for 
the degree of quantum coherence in the ground state of 
the system. Similar to the ferromagnetic case discussed 
in I, the magnetization current circulating a Heisenberg 
antiferromagnet generates an electric dipole field. Due 
to screening effects, the corresponding voltage drop may 
be rather difficult to detect experimentally. However, in 
view of the rapid development of the field of spintronics, 
it does not seem unreasonable to expect that before the 
end of this decade new experimental techniques will be 
available to detect persistent spin currents in Heisenberg 
rings. 

We thank Florian Meier and Bernd Wolf for discus- 
sions. This work was supported by the DFG via Forscher- 
gruppe FOR 412, Project No. KO 1442/5-1. 

APPENDIX A: CROWN-SHAPED MAGNETIC 
FIELD 

Here we determine the classical ground state of an an- 
tiferromagnetic ring in the crown-shaped magnetic field 
depicted in Fig.^ i.e., a field with 

hi = ft,[sin('i9) cos(27rZj/L), sin(i9) sin(27r/i/L), cos(i?)] 

(Al) 




FIG. 4; Definition of angles for the description of the classi- 
cal ground state of an antiferromagnetic Heisenberg ring in a 
crown-shaped magnetic field. 

For very strong magnetic fields the classical unit vectors 
rrii will be aligned parallel to the field and the ground 
state will have the full rotational symmetry of the applied 
field. Below a critical spin-flip field hc{'9) it will be ener- 
getically favorable to form two sublattices with different 
angles to the z-axis (see Fig. 0)). Introducing the 

relative and average angles 

5=\{Kn-^Z) (A2) 
a minimum of the classical energy is reached for 

sin(i? - I?) cos((5) = -e_sin(2i?) (A3) 

h 

JQ 

cos(i5 - ?9) sin((5) = +—e+sm{2S), (A4) 

where we have defined e± = 1 ± cos(27r/A^). For very 
strong magnetic fields h > he, we have 6 = and 
Eq. (|A3() reduces to its ferromagnetic analogue (see Eq. 
(14) in I). For 6^0 the two equations can be combined 
to give 

sm{2{§ -§)) = -( — ] sm\2n/N) sin(2^) . (A5) 

Thus for large rings the magnetic field h ^ JS/N nec- 
essary to produce an inhomogeneous classical ground 
state is well below the spin-flip field he ~ JS. For 
h JS/N ^ JS we have S ~ 7r/2 and the classical 
ground state locally resembles a Neel state as assumed in 
SecHH 



APPENDIX B: CALCULATION OF SPIN 
CURRENT 

In this appendix we present some details of the calcu- 
lation of the functional form of the spin current at zero 
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temperature. Eq. (|23|l can be written in the form 

where f{k) — d/d{ka){ek/2JS) is a periodic function of 
its argument with f{x + tt) = f{x) = —f{—x). Conse- 
quently Im is a periodic function of and we can proceed 
by calculating its Fourier coefficients Ci, in 

By appropriate substitution in the corresponding integral 
one eliminates the finite momentum sum and obtains 

a = -^<?.^.(-A-2) , (B3) 



where .92/+! is zero and 52/ is given by 



92i{z) = 




Using an asymptotic expansion of the hypergeometric 
function^"* for large /, we derive an expression for the spin 
current for large system sizes. In particular, the condi- 
tion that the leading term in the expansion be sufficient 
can be shown to be equivalent to L/^ C'(l)- After 
some algebra we obtain Eq. (|29|l . 
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